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GeV Emission from Collisional Magnetized Gamma Ray Bursts 

P. Meszaros 1 and M.J. Rees 2 

ABSTRACT 

Magnetic fields may play a dominant role in gamma-ray bursts, and recent 
observations by Fermi indicate that GeV radiation, when detected, arrives de- 
layed by seconds from the onset of the MeV component. Motivated by this, we 
discuss a magnetically dominated jet model where both magnetic dissipation and 
nuclear collisions are important. We show that, for parameters typical of the ob- 
served bursts, such a model involving a realistic jet structure can reproduce the 
general features of the MeV and a separate GeV radiation component, including 
the time delay between the two. The model also predicts a multi-GeV neutrino 
component. 

Subject headings: gamma-rays: bursts — cosmology — stars: population III - 
jets: magnetized - radiation mechanisms: non-thermal 



1. Introduction 

Recent observations with the Fermi satellite have shown the presence of GeV emission 
in GRB bursts, delayed by seconds relative to the usual MeV radiation, e.g. Abdo et al. 
(2009a); Ackermann et al. (2010). In some bursts, this delayed GeV radiation appears as a 
distinct spectral component (Ackermann et al. 2010; Abdo et al. 2009b). There have also 
been indications that the outflow in such bursts may be magnetized, e.g. Zhang and Pe'er 
(2009); Fan (2009). Also on the theoretical side, there has been an increasing interest in 
magnetically dominated but baryon loaded GRB models, e.g. Thompson (1994); Drenkhahn 
and Spruit (2002); Lyutikov and Blandford (2003); Metzger et al. (2011); McKinney and 
Uzdensky (2011). 
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The fast rotating central engine, whether a black hole or a magnetar, can lead to a highly 
magnetized outflow which is initially Poynting dominated, with a sub-dominant baryon load 
at the base of the flow r such that 77 = (£em/ Pb c2 ) ^ 1> where is electromagnetic 
(mainly magnetic) energy density, and p& is baryon mass density. If the magnetic field in 
the jet is striped, i.e. it has alternating polarity, either as a result of flow instabilities in a 
collapsar outflow or in oblique magnetar outflows, reconnection and dissipation can occur 
both below and above the photosphere, as discussed by Giannios and Spruit (2007); Metzger 
et al. (2011); McKinney and Uzdensky (2011), which can give rise to efficient radiation in the 
MeV band. Dissipative photospheres have also been proposed in magnetically sub-dominant 
baryon loaded outflows, as the source of the MeV radiation, e.g. Rees and Meszaros (2005); 
Pe'er et al. (2006); Ryde et al. (2010); Pe'er et al (2011); Beloborodov (2011), while the role 
of nuclear collisions in producing dissipation, high energy photons and neutrinos has been 
discussed by, e.g. Meszaros and Rees (2000a); Rossi et al. (2006); Ioka (2010); Beloborodov 
(2010) and others. 

In this paper, we analyze a magnetically dominated and baryon loaded outflow, and 
discuss the combined effects of dissipation from magnetic reconnection and nuclear collisions 
in the context of the magnetically dominated outflow dynamics, which differs from the usual 
hydrodynamic regime dynamics. We show that such outflows can produce the usual MeV 
spectrum, as well as a GeV component which arrives with a delay of seconds after the MeV 
component, for a minimally (realistic) structured jet configuration. 



Dynamics - If the magnetic field is striped and the jet is one-dimensional 3 the outflow accel- 
erates (Drenkhahn 2002; Granot et al. 2011; Metzger et al. 2011) roughly as 



where r = 10 7 r j is the base of the outflow and r sat ~ r r/ 3 is the radius at which the 
bulk Lorentz factor T has saturated, the outflow having become approximately matter dom- 
inated. One way to roughly understand this dynamic behavior is as follows. The comoving 
reconnection time is t' r ~ A'/i^ = Y\/v' r ~ Teir /e2VA oc T, where A = eir is the approxi- 
mate radial lengthscale of field polarity changes, e± ~ 1, v' r ~ tiv' A is reconnection velocity, 



2. Magnetic Dynamics, Nuclear Collisions and Photosphere 




for r < r sat 
for r > r sat 



(1) 



3 For more general magnetic jet geometries and acceleration laws see text below eq.(13) and Komissarov 
et al. (2009); Tchekhovskoy et al. (2010); McKinney and Uzdensky (2011). 
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v ' A ~ c is the Alfven velocity, with, e.g., e2 ^ 0.1 and comoving (lab) frame quantities are 
denoted as primed (unprimed). The dynamic expansion time is t' ex ~ r/cT. The comoving 
internal energy per baryon of the plasma (which is mainly magnetic, in the acceleration 
phase) can be expressed as 7'. The outflow will be self-similar, and 7' > 1, a dimensionless 
quantity, can be expressed as a dimensionless combination of dimensional quantities in the 
flow. There are only two physically independent dynamical quantities, t' r and t' ex , suggesting 
7' oc t' r /t' ex oc T 2 /r. The conversion of internal energy during the expansion is what leads to 
the growth of the bulk kinetic energy, and energy conservation requires 7T ~ const, oc T 3 /r, 
or T oc r 1 / 3 . The final bulk Lorentz factor is the same r ~ rj as in the pure hydrodynamic 
regime, but the growth is slower, because at each e-folding in radius only part of the magnetic 
internal energy gets converted directly into kinetic energy, the other part gets converted first 
into thermal energy, and only after another e-folding into kinetic energy. 

Longitudinal pion optical depth.- In a simple flow which initially moves with a single bulk 
Lorentz factor, nuclear inelastic collisions couple n and p with a cross section a nuc ~ 
a n (c/v re i), where v re i is the relative drift velocity between n and p and (T T ~3x 1CT 26 cm 2 . 
When v re i — > c, the nuclear collisions become inelastic resulting in the production of pions, 
muons, and eventually pairs and and neutrinos. The latter escape, but ~ 1/4 of the 
collision energy ends up in pairs and photons, which are trapped in the optically thick flow, 
adding to its thermal energy. 

In such simple homogeneous flows, the nuclear collisions become inelastic (resulting in 
pions) when the comoving expansion time t' ex ~ r/cT becomes shorter than the comoving 
collision time t' nuc ~ l/(n p a n c), or 7^ y ~ n' p a nuc r /T ~ 1. Here n' p is the comoving proton 
density 

where L is the total kinetic luminosity, and x = n p jn\ ) = n p j (n p + n n ) is the proton fraction 
of the baryon density n& = n p + n n . The radius at which pion production occurs is 

Lo~ 7T x 1 g x 



where 



r Airm p c 3 r T]T 2 n rjT 2 



ti-T^, (3) 



(4) 



and we assume parameters comparable to those deduced for typical bright LAT bursts, e.g. 
GRB 080916C (Abdo et al. 2009a), L ~ 10 54 L 54 , V ^ 600^, 9 ~ 10" 2 . 



The nuclear collisions can become inelastic in the regime T ~ (r /tq) 1 ^ 3 < rj, in which 
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case from eqs. (1, 3) the longitudinal pionization radius is 

3/5 



^) for rv<r sat , i.e. rj > ^ = O.ftr&V 
V J 



(5) 



The pionization radius and r]^ of eqs. (5,4) scale with L, r differently than in the 
corresponding hydrodynamic outflows, e.g. Bahcall and Meszaros (2000) and others, due to 
using here the magnetized outflow dynamics (1) instead of T oc r. Beyond r n , the neutrons 
decouple from the protons and start to coast with a terminal bulk Lorentz factor 

T n ~ x^tjMtj) 1 / 5 ~ 80 l}£rj£l* x x J*, (6) 

while the protons continue to accelerate as Y oc r 1 / 3 until reaching the coasting value r ~ r\. 
If on the other hand the flow Lorentz factor saturates before the protons and neutrons 
decouple, i.e. before the nuclear collisions become inelastic, which occurs if r\ < r^x 1 / 6 , the 
protons and neutrons continue to coast together at T ~ r\ = constant and pionization never 
occurs. 

Photon optical depth- If the only electrons in the outflow are those associated with the 
protons, the Thomson photosphere occurs at 

r P h - ( LaT ) 1 - ^ 1 (7) 
r \4irm p c 3 r J r/T 2 r/T 2 



where 



/ T \ 1/6 / \ 1/6 

For rj > t)t (rj < Tfr) the photosphere occurs at r < r sat (r > r sat ), at a radius 



r_ph = f rf T {j] T /rif/ h for r < r sat (77 > r? r ); 
r o I Vt(Vt/v) 3 for r > r sat (?7 < ?7t)- 



(9) 



The scaling, as well as the definition of r/ T again differs from those in e.g. Meszaros and Rees 
(2000b) or Beloborodov (2010), due to our using here the magnetized outflow dynamics. 
Thus, the photosphere is typically larger then the longitudinal pionosphere by a factor 

^~f^V /5 ~6.4. (10) 



The inelastic nuclear collisions produce additional leptons which can contribute to the 
Thompson opacity. For a strongly magnetized outflow the initial secondary leptons would 
cool very fast and the cascade efficiency drops as the pair production approaches its threshold, 
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where in the limit (Beloborodov 2010; Vurm et al. 2011) each n,p collision adds at most one 
further lepton. Since T n oc r~ 5 / 3 drops rapidly, most neutrons will suffer only one inelastic 
collision, increasing the electron density by a factor ~ 2. In addition, the reconnection 
process can also produce pairs, whose effects could be significant (McKinney and Uzdensky 
2011; Metzger et al. 2011). Assuming that the two effects contribute a comparable change 
in the lepton density, this increases L by a factor > 4 in eqs.(7,8, 9). Thus, the photospheric 
radius could be a factor > 4 3 / 5 ~ 2.3 larger than its value without pairs in eqs. (7,9), i.e. 
r P h > 10 - 157V 

Transverse pion optical depth - Neutrons, if present outside the jet, can drift sideways into 
the jet channel unhindered by magnetic fields. The lab-frame transverse nuclear collision 
optical depth of a jet of opening half-angle 9 at radius r is 

r n ,± ~ n p a n r9 = — (11) 

\47rm p c d ro / r r] 

and the jet becomes transversely optically thin at 

^=,°*. (12) 
r r] 

Thus, for a jet opening half-angle 6 = 10~ 2 #_2, we have the following radii, 

r ml2r3/5 -3/5 3/5 -3/5 

r n ~ 6 x 10 U L 5 ' 4 r y x [ 5 n 60( j cm, 

r ph ~ 6 x lO^Lgfr^ 75 ^ 3 / 5 cm, 

r*,± ^ 9 x 10 14 L 54 r7 60 1 6'_2 cm, 

r sat ~ 2 x 10 15 r , 7 ^Ioo cm . (13) 

Transverse jet structure and neutron drift- In a realistic jet the properties vary also in 
the transverse direction. Hydrodynamical simulations (e.g. (Morsony et al. 2010) indi- 
cate a Lorentz factor tapering off towards the edges, while MHD force-free simulations 
(Tchekhovskoy et al. 2008) show models where, depending on the stellar pressure profile 
and the magnetic symmetries, the jet opening angle can initially narrow with increasing 
radius, with an inner jet core accelerating initially more slowly than the outer jet portions, 
but before exiting the star the inner jet becomes faster than the outer sheath, after which 
the jet opening angle remains constant. As a simple idealized model, we can consider the 
transverse structure outside the star as a two-step jet, consisting of an inner jet core with, 
e.g., a nominal constant half-angle 6 = 10~ 2 #_ 2 with r\ = 600i] 600 , and a slower 4 outer jet or 



4 Depending on various parameters an alternative model might have outside the star a slower inner jet 
and a faster sheath, which could lead to results similar to those discussed below; we do not consider this 
case here. 
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sheath with, say, r] out ~ 100, where both inner core and outer sheath have been populated 
with protons and free neutrons already near the black hole, where most nuclei get photo- 
dissociated (Beloborodov 2003; Metzger et al. 2008). One can expect neutrons from the 
outer sheath to drift into the jet core, and for i] out & 10 2 ~ 6~ l the neutrons can penetrate 
transversely the entire jet core. Also, the relative radial Lorentz factor ratio between indrift 
neutrons and jet core baryons (either p or n) is T re i ~ r)/r) out > 1, so the collisions will be 
inelastic, leading again to pions. 

Over a timescale t the thermal neutrons from an outer sheath where the lab-frame 
neutron density is n n and the thermal neutron velocity is Vt will diffuse transversely into the 
jet, their flux being 

'v t l\ 1/2 n ( n n c\ 1/2 



*-~{f) n "~ p '{tt) ' (14) 

where £ ~ 1/ (n n a nuc ) , a nuc ~ a n (c/v t ) and (3 t = v t /c. Taking a jet transverse area of 
A ~ n9r 2 , on a timescale t the total number of neutrons drifting transversely into the jet 
core will be, roughly, 

N n , ± ~ vrflr 2 n t ~ ir6r 2 (3 t {n n ct/ a^) 1 ' 2 , (15) 

while over the same timescale the number of baryons (n and p) passing longitudinally through 
the jet core is 

iV 6j || ~ 7[6r 2 n h ct ~ nQr 2 (L/47Tr 2 m p c 3 i])ct (16) 

The number of collisions is maximized when iV ni ^ > N b ^, that is f3 t (rinCt/a^) 1 ^ 2 > (ct/a^) (rj®r /r 2 ri), 
and when the jet transverse pion optical depth T nt ± ~ 1, i.e. r ~ r n ^±. For the first condi- 
tion, we can estimate the outer sheath neutron thermal velocity f3 t since we know that the 
comoving baryon temperature becomes trans-relativistic (T ~ GeV) at the saturation radius 
i"sat,out = T'ojVout = 10 13 cm ) where Vt ~ c, after which they coast with with T t oc r -2 / 3 , so at 
r-K,± we have f3 2 = (v t /c) 2 ~ {kT t /m p c 2 ) = {r sa t,out/'r-n^±)~' 2 ^ ~ 10 _4 ^ 3 These two conditions 
are then met when at r ~ r n ± the neutron density in the outer sheath is 



n 



n > ~ 4 x 10%?^^^ cm" 3 (17) 



Vl 2 r 2 e*(3-, 



over a timescale t. The density in eq. (17) is comparable to what might be expected in an 
outer jet sheath. 



3. Spectrum Formation 

GeV radiation- The 'prompt' component of the GeV radiation is expected to arise from 
the transverse drift nuclear collision mechanism at r\ j_ discussed above. The 77 photo- 



-7- 



sphere to E — 10 GeV photons can be roughly estimated as (Beloborodov 2010) r 7J (E) ~ 
10 15 (£/10 GeV) 3 / 2 r/ 6 - 5 L54 cm , so that 

r ph < r 77 (10 GeV) £ r W)± , (18) 

the right inequality holding for r] 600 < 1. Thus, at r^j., which is at or outside r 77 for r] 600 
multi-GeV photons will be copiously produced by transverse indrift neutron collisions with jet 
core baryons, and these photons can escape unhindered. The energy in this GeV component, 
for the conditions discussed in the previous section, is a significant fraction of the kinetic 
luminosity. Most pions will be formed approximately at rest in the frame of the jet, and the 
processes p,N—> 7r + ,7r~,7r°, where N = (p,n), occur in approximately equal ratios. Thus 
1 /3 of the collisions on average yield two photons from tc° — > 27 decays at an observer frame 
energy which is broadly distributed around a central energy e^' 1 - ~ r)(m n c 2 /2)/(l + z) ~ 
7 r76oo (5/1 + z) GeV, extending up to twice that value and down to a fraction of it. 

GeV time delay - The relative time delay between between MeV photons from the central 
photosphere (see below) and the observed onset of the GeV emission, which arises at r n ,± 
from neutrons from an outer jet sheath moving with bulk Lorentz factor r] out which drift into 
the inner jet at the radius r n> ± where their effect becomes significant is 

At ~ ~ 1 L ^ r ojVmoVoiioo s > ( 19 ) 

zc 'lout 

in the lab frame at the source. In the observer frame this is lengthened by a factor (1 + z). 
For a burst at z ~ 4 such as GRB 0809 16C this is similar to the observed GeV-MeV delay 
of ~ 5 seconds observed by Fermi in that object. 

MeV radiation- The MeV (Band) spectrum may be expected to arise from the dissipa- 
tive photosphere of the inner jet at r ph ~ 6 x 10 13 cm (eq.[13]). This MeV spectrum may 
be attributed to synchrotron radiation following magnetic or shock dissipation near r ph , 
with typical bulk Lorentz factors T r ~ 1. This would result in random magnetic fields 
B' ~ (327re B m p c 2 n^) 1 / 2 r r , where n' b = Lj {A-Kr^ h riY ph ) is the baryon comoving density at r ph , 

where for the nominal jet parameters used we have T p h = {rph/ro) 1 ^ 3 ~ 180L5{ 5 r 7 / ' 15 77 60 q / ' 5 , 
so n' b ~ 5 x 10 12 L 54 2 ^Vo 6 / 15 %oo and B' ~ 4 x 10 5 L b l' ' 5 >15 'ilH^e^T r G. The same mag- 
netic dissipation leading to semi-relativistic shocks with comoving bulk Lorentz factors 
T r ~ 1 near r p h will accelerate electrons to average non-thermal random Lorentz factors of 
7g ~ {m p /m e )Y r ~ 10 3 r r , assuming equipartition between baryons and electrons, leading to 
an observer frame synchrotron radiation peak at e p S y = e' sy T ph (l + z)^ 1 ~ 0.3 r l Jj e l ^ 2 Tl(h / 1 + 
z) MeV. This non-thermal component will dominate over a much weaker photospheric ther- 
mal component at kT t h ~ 1 keV. Comptonization may boost a fraction of the synchrotron 
photons to higher energies, but using the expression above for r 77 a cutoff is obtained at 
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e 77 ~ 3 MeVL^ r ^ 180 . Thus, the photosphere allows escape of MeV but not GeV pho- 
tons. A decreasing fraction of increasing energy photons could arise from dissipation at 
larger radii (Beloborodov 2010), as r 77 increases, but significant GeV escape compatible 
with a second spectral component is only expected from radii where transverse neutron drift 
is important. 

Outer jet photosphere - The outer jet has a photosphere beyond its saturation radius, at 
r P h,out ~ 2.5 x 10 14 L 54 ?7^ u 3 tiloo cm. If there was dissipation at this outer photosphere, a calcula- 
tion similar to the inner jet's would give a synchrotron peak at £ S y,out ~ 25 L 54 1 ^ 2 e]/ 2 r 3 7/^ 100 (5/l+ 
z) keV. However, magnetic dissipation is not expected to occur above the saturation radius, 
and neither does nuclear collisional dissipation, since the n,p components do not decouple 
for this r) < rj^ ~ 130. Thus, unless some other form of dissipation occurs near the outer pho- 
tosphere, it would have a quasi-blackbody spectrum with kT ~ 0.2 L 5 ^ 2 rlji]l^ 100 (5/1 + 
z) keV and a luminosity L pKout ~ e rad L (r sat /r ph ) 2/3 ^ 10~ 2 (e m(i /0.5)Lo, well below the 
kinetic luminosity L . 

Internal shocks- In principle these would not be expected if the magnetization remains 
significant beyond the saturation radius (c.f. Zhang and Yan (2011)). If they did occur, the 
earliest would need to be at r IS ~ 2ct v rf ^2x lO 15 ^-^^ cm ^ r sat , where t v = lO -1 ^,-! s 
is the variability time, and the usual arguments would lead to a synchrotron peak at e 1 ^ ^ 
8 Ll{ 2 efeX'-iVmo^r(5/l + z) keV. 

External shocks and afterglow.- Other radiation components are due to the external shocks 
of the inner and outer jets. Forward shocks (FS) are guaranteed, while reverse shocks (RS) 
are less certain, depending on the residual magnetization a value at the external shock. For 
an outflow (prompt burst) duration t = 10t\ s in the burst frame, in an external medium 
of density no cm -3 , the inner jet external shock occurs at a deceleration radius r d ~ 1.6 x 
10 17 L5{ 3 t^ 3 n ~ 1/ ' 3 ?7 60 2 / 3 cm at an observer time t d ~ hQL\^t l i Z n^ 1 ^ r]^ 3 (\l + z\/h) s. With 
the usual assumption about magnetic field and electron energy equipartition, and ~ rj/2 
at r d , the FS synchrotron spectrum of the inner jet peaks at efy ~ 2e] 3 / ' 2 eg?7g 00 (5/l-|-z) MeV. 
For the usual power law electron distribution accelerated at the shock, this synchrotron 
spectrum extends as a power law into the GeV range (Ghisellini et al. 2010; Kumar and 
Barniol 2010), providing the long term GeV afterglow, the external shock being well outside 
r 77 (10 GeV). The earliest prompt GeV radiation, however, appears to require a separate 
origin, as considered by Toma et al. (2011); He et al. (2011); De Pasquale et al. (2010); Corsi 
et al. (2010) in the context of non- magnetic outflows. 

The outer jet deceleration radius is r d)OUt ~ 5.3 x L^t^n^ 1 ^ 3 ^ioo i0 ut cm & t an ob- 
server time t dfi ut ~ 2 x 10 3 Ll{ 3 t^ 3 n ^ 3 ?7]~oo!Lt([l +z]/5) s, and the outer jet FS synchrotron 
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spectral peak is at ef y s out ~ 2 e B 7 e^ t>100 (5/l + 2) keV. 

If the magnetization parameter of the ejecta a has become small at r d) a reverse shock 
may form (Mimica et al. 2009; Zhang and Yan 2011). With the usual assumptions for an 
equipartition magnetic field in a reverse shock with T r ~ 1 and 7g ~ e e 10 3 r r the inner jet RS 
synchrotron peak occurs at ~ 70 L 1 J^t l ^nJ 2 e^ 2 e 2 e r^r]l^ > {b/l + z) eV, while the outer 
jet reverse shock synchrotron spectrum peaks at s^ out ~ 20 eV, with the same scaling. 

Neutrino spectrum - Also, as a result of the charged pion decays, 1/3 of the collisions yield 
on average a muon neutrino of energy broadly centered around e v ~ 7 r] 60 o(5/l + z) GeV 
(again extending up to twice that value and down to a fraction of it), and a muon neutrino 
as well as an electron neutrino centered at energy e u ~ 3 r/ 6 oo(5/l + z) GeV, while another 
1/3 of the collisions leads to the corresponding antineutrinos of similar energies. These are 
in the sensitivity range of the Deep Core array of the IceCube neutrino detector. 

4. Discussion 

We have presented a GRB model based on a magnetically dominated, baryon loaded 
outflow, where both magnetic reconnection and nuclear p, n collisions lead to dissipation 
around and above a photon photosphere. We have discussed the dynamics of the magne- 
tized jets based on a simplified reconnection prescription and an idealized jet transverse 
structure, deriving the radii at which the photon scattering photosphere occurs as well as 
the radii at which radial and transverse drifts between proton and neutron components lead 
to nuclear collisional dissipation, in the context of the magnetic jet dynamics. Although 
the detailed structure and dynamics could depend on various assumptions about the stellar 
envelope pressure profile, the magnetic field configuration, etc. (Tchekhovskoy et al. 2008; 
McKinney and Uzdensky 2011; Metzger et al. 2011), the present model qualitatively captures 
a number of the essential features. In particular, the characteristic radii and their scalings 
differ from those in a purely hydrodynamical jet. Both in this and in more general magnetic 
models a slower acceleration rate generically leads to larger photospheric and dissipation 
radii relative to the hydrodynamic case, in line with previous hints from Fermi observational 
results indicating larger emission radii. 

Using parameters inferred from Fermi observations of bursts showing significant emission 
in the LAT instrument in the GeV energy range, which typically are of very high luminosity 
and higher than average Lorentz factor (Zhang et al. 2011), we discussed the gross properties 
of the MeV and GeV photon production spectra. For a simplified structured jet model 
consisting of a central faster jet core and a slower outer jet surrounding it, we showed that for 
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reasonable parameters a significant GeV component is produced, with a time delay of seconds 
relative to the onset of the MeV component, in rough agreement with the observations of 
GRB 080916C, e.g. (Abdo et al. 2009a). Related to the photon emission, a multi-GeV 
neutrino emission component is also expected, whose luminosity is comparable to that of 
the multi-GeV photon emission. More detailed calculations, which are beyond the object of 
this letter, will be needed in order to make detailed spectral fits to specific bursts. 

We thank NASA NNX09AL40G, NSF PHY-0757155, the Royal Society, and the Insti- 
tute of Astronomy, Cambridge, for partial support, and Dr. K. Toma and the referee for 
useful comments. 
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